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ABSTRACT 

The N2B hybrid wing body aircraft was conceptually designed to meet environmental and performance goals for the 
N+2 generation transport set by the subsonic fixed wing project. In this study, flow fields around the N2B configuration 
is simulated using a Reynolds-averaged Navier-Stokes flow solver using unstructured meshes. Boundary conditions at 
engine fan face and nozzle exhaust planes are provided by response surfaces of the NPSS thermodynamic engine cycle 
model. The present flow simulations reveal challenging design issues arising from boundary layer ingestion offset inlet 
and nacelle-airframe interference. The N2B configuration can be a good test bed for application of multidisciplinary 
design optimization technology. 


Nomenclature 

ADP = Aerodynamic Design Point 

AIP = Aerodynamic Interface Plane 

BLI = Boundary Layer Ingestion 

BPR = Bypass Ratio 

HWB = Hybrid Wing Body Aircraft 

ISA = International Standard Atmosphere 

LTO ^Landing and take off 

NPSS = Numerical Propulsion System Simulation 

OML = Outer Mold Line 

OPR = Overall pressure ratio 

RTO = Rolling Take off 

SFC = Specific fuel consumption 

SFW = Subsonic Fixed Wing 

SLS = Sea Level Static 

77 i = Mass flow rate 

P t = Total pressure 

T t = Total temperature 

U = Streamwise velocity magnitude 

U e = Streamwise velocity magnitude at a boundary layer edge 

77 = mean total pressure recovery at fan face 
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I. Introduction 


I ^he Subsonic Fixed Wing (SFW) project of NAS A Fundamental Aeronautics Program has focused on research of 
next generations of civil aircraft. For each time frame of near-, mid- and far-term periods, specific reduction goals 
in noise, emission, and fuel bum have been defined from the current state of the art of aviation technology. [1] (See 
Table 1), To achieve the aggressive goals especially for the N+2 and N+3 time frames, innovative design concepts in 
airframe and propulsion need to be explored. 


Table 1. NASA subsonic fixed wing project goals.[l] 


TECHNOLOGY 

TECHNOLOGY GENERATIONS 
(Technology Readiness Level = 4-6) 

BENEFITS* 

N+1 (2015) 

N+2 (2020**) 

N+3 (2025) 

Noise 

(cum margin rel. to Stage 4) 

-32 dB 

-42 dB 

-71 dB 

LTO NOx Emissions 
(rel. to CAEP 6) 

-60% 

-75% 

-80% 

Cruise NOx Emissions 
(rel. to 2005 best in class) 

-55% 

-70% 

-80% 

Aircraft Fuel/Energy Consumption* 
(rel. to 2005 best in class) 

-33% 

-50% 

-60% 


* Projected benefits once technologies are matured and implemented by industry. Benefits vary by vehicle size and mission. N+1 and N+3 values 
are referenced to a 737-800 with CFM56-7B engines, N+2 values are referenced to a 777-200 with GE90 engines 
** ERA'S time-phased approach includes advancing "long-pole" technologies to TRL 6 by 2015 
$ C0 2 emission benefits dependent on life-cycle C0 2e per MJ for fuel and/or energy source used 


The hybrid wing body (HWB) aircraft, although not new, is an alternative design concept to the conventional 
tube-and-wing aircraft. [2] It has the fuselage and wings integrated into a flying wing, which has better aerodynamic 
efficiency than the tube-and-wing type configuration. Also, the wide airframe body of HWB configuration is beneficial 
for shielding downward-propagating noise from engines installed above the aircraft. 

The Boeing Company studied HWB configurations meeting NASA SFW N+2 goals through a contract 
sponsored by NASA’s Subsonic Fixed Wing (SFW) project of the Fundamental Aeronautics Program. [3] The “Silent 
Aircraft” SAX-40 [4, 5], a non-proprietary configuration was chosen as a starting point. Two HWB configurations 
were considered in the study; one is the N2A model employing conventional podded engines and the other is the N2B 
model, which uses embedded engines and boundary layer ingestion (BLI) offset inlets. The SAX-40, N2A and N2B 
configurations are illustrated in Fig.l. Figure 2 shows the embedded engine having three fans and the variable area 
thrust vectoring/reversing nozzle concept adopted in N2B. Compared to conventional pylon-mounted engines, 
embedded engines with BLI offset inlets allow reduced ram drag, wetted area, structural weight, and noise. 
Disadvantages of the embedded engines are higher flow distortion and reduced pressure recovery at engine faces due to 
the BLI offset inlet. 

Owing to its nature of tightly coupled disciplines, the N2B configuration offers an excellent test bed for 
conducting multidisciplinary analysis/design optimization, e.g., with inclusion of aerodynamics, propulsion, and 
acoustics. High fidelity CFD based analysis/design is essential for accurate performance evaluation/improvement of 
propulsion-airframe integration. In the present study, as a first step toward MDAO of the HWB configuration, we 
conducted flow simulations of the N2B configuration to understand effects of the propulsion-airframe integration on 
the flow field around the vehicle and the embedded engine. 

The remainder of this paper is organized as follows. First, the N2B configuration is reviewed and response 
surfaces for thermodynamic cycle models of the embedded engine are generated. In Section 3, numerical approaches 
for the flow simulation are explained including the flow solver, mesh generation and boundary conditions. Simulation 
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results for the N2B aircraft are illustrated in Section 4. Summary and conclusions are followed in Section 5. Finally, a 
validation study for the flow solver on a wind tunnel test case of a BLI offset inlet is presented in Appendix. 




(a) SAX-40 



(b) N2 A with podded engines 



(c) N2B with embedded engines 




Fig. 1 HWB aircraft configurations 
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(a) Tri-fan + core engine 



(b) Core engine + variable area thrust vectoring nozzle 
Fig. 2. Multiple fan embedded turbofan engine for N2B. [3] 


II. N2B Hybrid Wing Body Aircraft 


A. N2B design parameters [3] 

The N2B was suggested by Boeing as a cargo freighter aircraft and conceptually designed to have a 477,400 lb m 
MTOGW, 103,000 lb m payload and a 6,000 nm range. The cruise Mach number is 0.80 and initial cruise altitude is 
35,000 ft. Time to climb 31,000 ft is 0.29 hour. The span length is 213 ft, and the reference area is 9246 ft 2 . It has three 
engines, each of which has three fans and one core engine. 

B. Engine model 

Figure 3 shows the internal layout of the tri-fan embedded turbofan engine. The concept was based on the 
SAX-40 GRANTA engine. [5] Conceptual design of the embedded engine was conducted in NASA Glenn Research 
Center for 2020 technology level. [6] Each engine has a gas generator (core engine) which drives an inline fan and two 
additional outboard fans through a mechanical transmission system. The Aerodynamic Design Point (ADP) is Mach 
number 0.8 and height of 3 1 ,000ft. The required thrust per engine at ADP is 1 0,000 lb at ISA+0. The design fan 
pressure ratio is 1.50 at ADP. All the fans are of the same diameter and rotate at the same speed. The design bypass ratio 
(BPR) of the core engine is 3. 1 and overall effective BPR of the tri-fan engine is 1 1 .3 at ADP. The fan face Mach 
number is 0.674. The engine thermodynamic cycle was designed using the NPSS (Numerical Propulsion System 
Simulation) program [7, 8] in the study of Ref. [6]. Some details of N2B embedded engine cycle information are given 
in Table 2. 
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267 " 



Fig.3 Internal layout of the imbedded engine. [6] 


Table 2. NASAN2B embedded engine cycle information. [3] 



Net Thrust (lb) 

49060 

36000 

10000 

8286 

7458 

SFC (Ib/hr/lb) 

0.288 

0.398 

0.564 

0.553 

0.556 

BPR (core only) 

3.2 

3.2 

3.1 

3.1 

3.3 

BPR (effective) 

11.5 

11.7 

11.3 

11.3 

11.8 

OPR 

45 

43 

46 

46 

43 

N1 (%) 

98 

97 

100 

100 

96 

N2 (%) 

99 

99 

100 

100 

99 

T3 (R) 

1706 

1710 

1476 

1430 

1398 

T4 (R) 

3460 

3460 

3010 

2920 

2844 

T41 (R) 

3310 

3310 

2876 

2789 

2716 

T49 (R) 

2460 

2460 

2113 

2044 

1986 


Total pressure recovery values at engine fan faces at ADP are estimated by Boeing from the mass averaged total 
pressure of the inlet capture flow at the inlet highlight using RANS flow analysis results on a clean wing configuration 
without engine installation assuming additional 0.6% total pressure loss inside the S-shaped inlet duct. [3] The recovery 
values are shown in Fig. 4 for three fan faces of the outer inlet. The center fan face recovery is composed of the bypass 
recovery of 0.96705 and the core recovery of 1 .00 assuming that uncontaminated flow is fed into the core engine, while 
the boundary layer flow is ingested into the two side fans and bypass duct of the central passage. The engine 
thermodynamic cycle design was conducted based on the recovery values at ADP. 
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Fig. 4 Input values of pressure recovery for engine thermodynamic cycle design at ADP. 


In the present study for computational fluid dynamics simulations of the N2B configuration, the embedded 
turbofan engines are replaced with an NPSS engine model as illustrated in Fig.5. Total pressure recoveries at engine fan 
faces are calculated by CFD and provided to the NPSS model and engine boundary conditions required for CFD 
analysis are obtained from the NPSS model as shown in Fig.5(c). 

In the flow simulation, an engine fan face is treated as a subsonic exit boundary condition with specified target 
mass flow rates and a nozzle exhaust plane is treated as a subsonic inflow boundary condition with specified total 
conditions and a mass flow rate. 



(a) Section view of the embedded engine core path. 



(b) Replacement of the engine with NPSS model. 


/ 


Pressure recovery at AIP 

f 


Flow Solver 


NPSS model 

V 


AIP and engine exit 

k 



boundary conditions 


(c) Coupling of flow solver and NPSS model. 
Fig. 5. Engine boundary conditions by NPSS model. 
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In actual implementation, the coupling of the NPSS model with the flow solver can be made in a direct way by 
integrating the flow solver with the NPSS program via a script file for exchanging data of the two codes. An alternative 
choice is to build algebraic response surface models for the engine cycle and inserting the algebraic models into the 
flow solver. [9] In the present study, we employ the response surface approach for simplicity of the coupling. Accuracy 
of the coupling is not affected by the response surface approach because the behavior of the NPSS engine cycle model 
can be fitted very accurately by simple algebraic equations as will be shown later. 

Engine parameters required for the present flow simulations are mass flow rate m, total pressure P t , total 
temperature T t at the nozzle exit of each flow passage. Input variables to the NPSS model are total pressure recoveries 
r] at fan faces. There are three input variables; 7j lt r\ 2 , Vs, f° r passages 1 (center passage), 2 and 3 (two side passages), 
respectively. If all the variables are independent and coupled together, we need nine functions, each of which is a 
function of the three input variables rj lt r\ 2 , Vs'- 

~ fmSjlv r l2’ r h\ p ti = — fniHi’TJi’th) 

= fm2(jl\’ r h> r h\ P t2 = /p2 (^15 ^2 5 ^ 3)5 P t2 ~ It2 (/7l ? ^2 ’ ^3 ) (p 
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Fig. 6 Sensitivity of engine parameters with respect to changes in recovery at fan faces. 
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Sensitivities of engine parameters are calculated for one percent variation of the recovery at fan faces. Figure 6 
shows some results of the sensitivity study. All other sensitivity terms not shown in Fig.6 are all zero. The results of the 
sensitivity study can be summarized as follows: 

1 ) m and P t are affected significantly by the recovery of its own passage only. 

2) T t of the central passage is affected by the recovery of the central passage, but the amount of sensitivity is 
negligible. 

3) P t , and T t of the central passage are slightly affected by recoveries of side passages because the turbine work 
of the core engine to drive side fans is affected by mass flow rates of the side passages, which depend on the 
recovery values at the side fan faces. However, the amount of sensitivity is negligible. 

4) The recovery of the central passage has no effects on m, P t , or T t of side passages. 

Through the sensitivity study, the number of response surface models can be reduced from nine to three; 

(fm> fpi> fp 2 ,z)i an d each model is a function of only the recovery of its own passage. 

T n = const, 

^2 = fm (jh )’ ^2 = fp2,3 (jl 2 )’ ^2 = COTlSt 
/77 3 = fm ( 7 7s )’ ^3 = fp2,3 {^3 )? ^3 = COnSt 2 3 


Data points for response surface models can be generated by conducting off-design simulation of the NPSS 
engine model for different inlet recovery values. In Fig.7, the resulting engine parameters of interest for airframe-engine 
integration perspective are plotted, and least squares curve fitting is applied to obtain algebraic response surface models. 




Fig. 7 Response surfaces for the NPSS engine model at the maximum power condition with Mach 0.8, h=35 k ft. 
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C. Sectional area distribution in the S-duct 

The one-dimensional section area distribution inside the S-duct has a big impact on the performance of an offset 
inlet. Three sections of the N2B offset inlet shown in Fig. 8 correspond to (1) the inlet throat section, (2) the position just 
after the vertical splitter dividing core and side passages, and (3) the engine fan face. The section area increases in the 
first half of the inlet duct from section (1) to section (2) and decreases in the second half of the duct from section (2) to 
section (3), a trend following an inlet design concept in Ref. [10] to reduce flow separation in the S duct. In Ref. [10], 
also discussed is to apply pre-compression (in front of the inlet) so as to reduce required diffusion amount inside the 
inlet. 



(a) Section geometry of center passage 
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(b) Area distribution of duct and spinner 
Fig.8 Section area distribution inside the S-duct. 


In the original CAD geometry of the N2B aircraft, the spinner at the fan hub was not included. The spinner 
geometry was added into the model in the present study to be consistent with the area and average Mach number at the 
fan face determined through the NPSS engine thermodynamic cycle design. As evident in Fig. 8 (b), the final area 
distribution after subtracting spinner area from the duct area has nearly the same magnitude of areas at sections 1 and 3, 
which means the net diffusion made through the duct is very little. In the design of section area distribution of an S-duct, 
the fan hub or spinner area needs to be considered for a better diffusion performance of the inlet duct. 
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III. Methodologies for Flow Analysis 


A. Flow Solver 

GO -flow [9], a finite- volume unstructured-grid Navier-Stokes solver, is used in the current flow simulations. The 
compressible Reynolds-averaged Navier-Stokes (RANS) equations are discretized by the cell-vertex finite volume 
method. Control volumes are non-overlapping dual cells constructed around each node. Each edge connecting two 
nodes is associated with an area vector of the control surface, at which flow fluxes are computed. For a second order 
accuracy, a linear reconstruction of the primitive gas dynamic variables inside the control volume is used in conjunction 
with a limiter. The inviscid flux is computed using approximate Riemann solvers. Turbulence effects are considered by 
using Mentor’s two equations Shear Stress Transport (SST) model [11, 12]. For the time integration, the Fower-Upper 
Symmetric Gauss Seidel (FU-SGS) implicit method is adopted [13]. Parallel processing is made by domain 
decomposition of the computational mesh and Message-Passing Interface (MPI). A validation study of the flow solver 
on a BFI offset inlet configuration is presented in Appendix. 

B. Grid Generation 

A triangular surface mesh is generated directly on Stereolithography (STF) data, which can be generated from a 
CAD model. Generation of a hybrid volume mesh is made by using advancing front/layer methods adopted in 
MEGG3D [14, 15]. Prism layers are generated near viscous walls, tetrahedral cells in the remaining computational 
domain and pyramid cells in between when necessary. Figure 9 (a) displays a surface mesh for the symmetric N2B 
configuration. Figure 9 (b) shows a top view of duct system splitting into three branches. The surface mesh near 
nozzles is shown in Fig.9 (c). As mentioned earlier for Fig.2 (b), the N2B aircraft has variable area thrust vectoring 
nozzles, which inevitably contain moving parts contacting each other or having very small gaps in between. The 
complex geometry of movable devices was simplified into clean exhaust nozzles as shown in Fig.9(c). If the thrust 
vectoring function needs to be simulated, the simplified surface geometry can be deformed for that purpose. 

The number of computational mesh points for the N2B configuration is about 13.3 million, and first nodes off the 
viscous walls are clustered to the wall so that the maximum y+ value is less than 2. 

C. Boundary Conditions 

Freestream conditions are for a cruise flight at Mach number of 0.8 and altitude of 35,000 ft. Freestream 
boundary conditions for turbulence valriables are set as k^/u^ = 10 -6 , co^L/u^ = 5 following recommendations 
given in Ref. [12]. 

For the engine fan face boundary conditions, we extrapolate the density and velocity components from inside of 
computational domain and impose a static back pressure to match the target mass flow rate. For the engine nozzle 
exhaust planes, the total pressure and total temperature are set and a static pressure is taken from downstream in the 
computational domain for subsonic inflow boundary conditions. 


IV. Flow Analysis Results for N2B 

The integrated airframe-nacelle configuration shows complex flow features. As shown in Fig. 10, flow over the 
upper surface accelerates around the outer surface of cowl lip higher than Mach 1.15, which may cause shock induced 
boundary layer separation. Figure 11 shows flows are separated inside of S-ducts and nacelle upper surfaces. 
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(a) Overall Surface mesh for symmetric half model. 



(b) Top view of S-duct and spinners at fan faces. 



(c) Near the nozzles. 


Fig. 9 Surface mesh for the N2B configuration. 
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Fig. 10. Supersonic envelop with M > 1.15. 



(a) Iso-contours with u velocity < -0.01 



(b) Oil flows on viscous walls. 
Fig. 11. Separation patterns 
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Fig. 12. Comparison of upper surface pressure contours of N2B and clean wing HWB. 
(left: clean wing, right: N2B). 


The upper surface pressure contours fo the N2B and clean wing HWB which does not have any engine installed 
are compared in Fig. 12. The installation of embedded engine nacelles on the upper surface causes higher pressure in 
front of the inlets due to the pre-compression, which is expected to cause a lift deficit for N2B. 



Fig. 13 Comparison of lift coefficients between clean wing and N2B. 
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y=O.Oi 





Fig. 14 Comparison of section pressure distributions between clean wing and N2B. 
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Fig. 14 Comparison of section pressure distributions between clean wing and N2B. (cont’d) 



(a) Total pressure and Mach contours at symmetric plane. 



(b) Total pressure and Mach contours at center plane of outer inlet. 


Fig. 15. Total pressure contours at inlet center sections. 
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The boundary layer ingestion can be visualized by total pressure contours at the center section of each inlet, as 
shown in Fig. 15. The center inlet has a thicker boundary layer ingested into the S-duct than the side inlet, because of the 
longer distance form the airframe leading edge to the inlet entrance. Mach contours in Fig. 15 also show flow separation 
patterns on the upper surface of the nacelle. It is noted that there is a big difference in surface slopes in front of the offset 
inlet at the symmetric plane and the center plane of the outer inlet. At the symmetry plane, the upper surface before 
entering the inlet is well aligned with the outer surface of the nacelle. On the other hand, at the center plane of outer 
inlet, the nacelle cowl angle has a quite large relative incidence angle to the incoming flow along the upoper surface. 
This relative angle of attack of the cowl lip is causing the leading edge separation on the upper surface of the nacelle of 
the outer inlet. This is due to the fact that the x/c (relative longitudinal location to the section airfoil chord length) of the 
inlet entrance is abount 77% for the symmetric plane and 61% for the center plane of the outer inlet. The difference of 
local x/c means a difference in local slopes of the section airfoil surfaces. Since a HWB aircraft tends to have a large 
variation of the section chord length, cowl lip angles of embedded engine nacelles should be designed to align with 
local inflow directions. 



(a) Location of streamwise section cuts 



X=0.718 


X=0.777 





X=0.74 



Pt: 0.7 0.8 0.9 1 


(b) Total pressure normalized by freestream total pressure 
Fig.16 Total pressure contours inside inlets. 
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Total pressure contours in the inlets long the streamwise coordinate x are shown in Fig. 16. The difference in the 
thickness of ingested boundary layers is clearly seen. The side inlet has left-right asymmetry in the contours. To explain 
the asymmetry, streamlines into the engine faces are visualized in Fig. 17, which shows relative side angle of the inflow 
to the side inlet. 



Fig. 17. Streamlines into the engine faces and out of engine nozzle exhausts. 


N2B 



Fig. 18 Comparison of boundary layer edge lines for N2B and clean wing configurations at the symmetric plane 

Boundary layer profiles obtained from CFD simulations of a clean wing configuration has been widely used for 
the sizing of the inlet capture area for BLI offset inlets and estimation of pressure recovery at fan faces of embedded 
engines[16, 3] The boundary layer thickess of the N2B and clean wing configurations are compared in Fig. 18. The 
boundary layer thickness in the figure is defined as the contour line of 99% total pressure relative to the freestream total 
pressure. The boundary layer thickness is very similar between the two configurations on the upper and lower surfaces. 
However, as the flow over the upper surface of N2B approaches to the inlet, the boundary layer thicknessof N2B gets 
thicker than that of the clean wing because of the adverse pressure gradient due to the existence of the inlet and the 
pre-compression in front of the inlet. 

In Fig. 19, compared are recoveries of the side inlet obtained by clean wing simulation results and the present 
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N2B analysis result. The recoveries estimated by Boeing [3] with clean wing analyses and used as input values to the 
NPSS engine models in Fig.4 are shown again in Fig. 19 (a), where the bypass and core recoveries for the central 
passage are mass-weighted averaged for the design BPR 3.1. The recoveries by the clean wing results and present N2B 
analysis results are in excellent agreement with differences less than 0.3%. Clean wing results for recoveries of the 
center inlet are not available. For the present N2B simulations, the recoveries for the center inlet are 0.94153, 0.95842 
and 0.94153 for left, center and right passages, respectively, which are about a couple of percents lower than the side 
engine recoveries in Fig. 19 due to the thicker boundary layer for the center engine inlet. 

All the simulation results above are obtained with stationary spinner at the hub of the engine fan. Actually, the 
spinner is rotating at the same speed with the fan. The effects of spinning on the flow field inside the inlet is checked 
here by conducting the flow simulation with corresponding wall velocity boundary conditions on the spinner wall. 
Total pressure contours at fan faces are compared without and with rotation of the spinners in Fig. 20. As can be noted, 
the effect of spinning on the flow is negligible except for the spinner boundary layer. 

In Fig.21, Mach contours at the symmetric plane are shown around the nozzle flow, which is choked at the throat. 
A separation zone is seen under the nozzle with very low local Mach number. 



(a) Estimated recoveries by clean wing analyses results + 0.6% inlet diffuser loss 



(b) The present simulation results 

Fig. 19. Comparison of pressure recovery at AIP between clean wing and N2B. (front view) 
Differences are -0.20, 0.07 and -0.27% for left, center and right passages, respectively. 


No spinning 



Pt: 0.7 0.8 0.9 1 


Fig. 20. Comparison of total pressure contours at fan faces with or w/o spinning hub 
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Fig. 21. Nozzle exhaust flow 


V. Concluding Remarks 

High fidelity flow simulations were conducted for the N2B Hybrid Wing Body configuration with 
airframe-propulsion integration effects. Thermodynamic cycle of turbofan engines proposed for the N2B aircraft were 
represented by the NPSS engine models. The simulation results revealed complex flow characteristics resulting from 
the tightly integrated airframe-propulsion system, which are very hard to estimate a priori and are seen to be of 
considerable departure from that of a clean wing without embedded engines. 

The simulation results show that the pre-compression in front of the embedded engine inlet caused a lift deficit, 
which is very critical in aircraft performance. Also, strong shock waves and flow separations are occurring on the 
nacelle upper surface. Local inflow angle effects such as side angle or incidence angle to the inlet cowl should also be 
considered in a detailed shape design of the airframe-propulsion system. The present flow simulation results on the 
hybrid wing body configuration are providing information and knowledge not only on the N2B aircraft but also on 
other HWB configurations with embedded engine concepts. 

Future work will be focused on shape design optimization of the HWB configuration in two aspects. The first one 
is improving flow quality at fan faces by a redesign of BLI offset inlet including the pre-compression region in front of 
the inlet. Another work is reshaping airframe-nacelle integration to minimize shock waves and flow separations outside 
of the nacelle. These two aspects of the shape design are not independent problems, but closely coupled with each other. 
Therefore, application of optimal shape design is essential to achieving the above goals. 
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APPENDIX: Validation of the Flow Solver for Boundary Ingestion Offset Inlet Diffuser 


For validation of the flow solver on configurations including BLI offset inlets, flow simulations are conducted 
for a flush-mounted offset inlet shown in Fig.Al and compared to the experimental results obtained by Owens et al.[17] 
at NASA Langley’s 0.3-Meter Transonic Cryogenic Tunnel, in which the inlet was flush-mounted on the tunnel 
sidewall. 


side view 


front view 



(a) Side and front views 



(b) CAD surface representation 
Fig.Al BLI offset inlet configuration: Inlet A model [17] 

Computational mesh for the BLI offset inlet is shown in Fig.A2. Outer boundaries are composed of inflow, side, 
top and outlflow planes. The outer boundaries are set as freestream boundary conditions. The side, top and outflow 
planes respectively are 50 inches away from the inlet throat centroid. The viscous flat plate length ahead of the inlet is 
adjusted to match the experimental bondary layer thickness at the inlet. The total number of mesh points is about 10 
million, and first nodes off the viscous wall are clustered to the wall so that maximum y+ value is less than 2. 
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(a) Surface mesh for the computational domain 



(c) Zoomed view near the inlet entrance 


Fig.A2 Computational mesh 

To ensure the numerical simulations are conducted at the same flow condition as the wind tunnel testing of the 
BLI offset inlet, several parameters need to be matched; the ffeestream Mach number, the Reynolds number, the 
boundary layer profile at the inlet, and the inlet mass flow rate. As mentioned earlier, the height of boundary layer 
profile is matched by adjusting the flat plate length ahead of the inlet. The inlet mass flow rate is controlled by varying a 
static back pressure at the AIR The freemstream Mach number is 0.85 and the Reynolds number based on AIP diameter, 
Re D = 3.8x1 0 6 , which is consistent with the experimental flow conditions in Ref. [17], where the Mach number was held 
at a constant 0.85 upstream of the inlet. Freestream boundary conditions for turbulence valriables are set as 
koo/ulz = 10" 6 , cOqoL/Uoo = 5 following recommendations given in Ref. [12]. 

Because the flow field has adverse pressure gradient as the flat plate boundary layer flow approaches the inlet 
entrance, just matching the boundary layer height does not necessarily mean the same boundary layer is ingested into 
the inlet. The amount of adverse pressure gradient due to the local flow conditions in the test section of the wind tunnel 
would affect more the detailed shape of the boundary layer profile ahead of the inlet. Figure A3 shows velocity profiles 
by the present numerical simulations and experimental results at the boundary layer rake location. The computaional 
velocity profile in a solid line has slightly higher verlocity in the bottom part of the boundary layer, which is consistent 
with other CFD analysis results in Ref. [19]. Lee et al.[18] compared boundary layer velocity profiles for Re D = 

2.2 x 10 6 and 3.8 x 10 6 and found the lower Reynolds number, Re D = 2.2x1 0 6 , results in a closer velocity profile to the 
experimental profile than the higher Reynolds number Re D = 3.8x1 0 6 . Also, Re D = 2.2x1 0 6 results in better 
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comparisons for pressure distributions with experimental data than Re D = 3.8xl0 6 does. Thus, it is deemed that the 
extent of matching the experimental boundary layer profile affects accuracy of computation of the pressure 
distributions inside the duct. 

CFD simulations for the BLI offset inlet test case in literature often uses a lower Reynolds number than the 
experimental Reynolds number Re D = 3.8xl0 6 probably to match the boundary layer velocity profile more closely to 
the experimental data. (e.g. Re D = 2.2x1 0 6 in Ref. [19]) However, the use of a lower Reynolds number also results in 
significantly lower total pressure recoveries at AIR [18, 19] In this study, instead of changing the Reynolds number, the 
static back pressure of the outflow boundary plane is adjusted to match the computational and experimental velocity 
profiles in the boundary layer as closely as possible. In the case of back pressure adjustment, the side and top boundary 
surfaces are treated as inviscid walls. In Fig. A3, the graph in a dashed line shows the velocity profile obtained by the 
present study with adjustment of the back pressure, which is lower than the freestream pressure. Note that the outflow 
plane and AIP are different exit boundary planes. 



Fig.A3 Comparison of boundary layer profiles. 

Exp: M 7 =0.843, Re D =3.3xl0 6 , CFD: M,=0.850, Re D =3.8xl0 6 
CFD-1: present simulation results without back pressure adjustment 
CFD-2: present simulation results with back pressure adjustment 


Figure A4 compares inlet centerline pressures from experimental and the present CFD results. The CFD results 
include two sets of data without and with the back pressure adjustment for boundary layer profile matching. The 
pressure distribution without the back pressure correction deviates upward from the experimental data. Meanwhile the 
pressure with the back pressure adjustment shows a very good comparison both on upper and lower surfaces inside the 
inlet. 
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x, inches (Tunnel Coordinates) 


Fig.A4 Comparison of inlet centerline pressure on the top and bottom of the BLI inlet. 

M^O.85, Re D = 3.8xl0 6 , Experimental data: wall correction 
(See Fig. A3 for definitions of CFD-1 and CFD-2) 

The inlet mass flow rate was represented by the area ratio Aq/A c where Ao is the area in the freestream flow that 
corresponds to the mass flow rate fed into the inlet: 

Po^o^o = P C U C A C (Al) 

The inlet capture area, A c , is defined by the inlet section area at the cowl highlight. As for an indicator for flow 
distortion at AIP, the SAE circumferential distortion descriptor, DPCP avg [21] is used to compare with results of 
experimental and other numerical studies in the literature. 

The inlet performance parameters obtained by the present flow simulations are compared with the experimental 
data in Table Al. The present results show good comparisons with experimental and other CFD results. Finally, 
contours of the local total pressure normalized by the freestream value are compared for experimental and other CFD 
results. The present results are calculating the location of the minimum total pressure region slightly higher than other 
results. Other than that, the overall trend is well simulated by the present CFD code. 


Table Al Comparison of performance parameters of the BLI inlet 
(See Fig. A3 for definitions of CFD-1 and CFD-2) 



A 0 /A c 

Distortion (DPCP avg ) 

Pressure recovery 

Reynolds number Re D 

Experiment [17] 

0.534 

0.054 

0.952 

3.8xl0 6 

Present CFD-1 

0.536 

0.052 

0.948 

3.8xl0 6 

Present CFD-2 

0.532 

0.054 

0.941 

3.8xl0 6 

CFD [18] 

0.533 

0.060 

0.956 

3.8xl0 6 

CFD [18] 

0.527 

0.063 

0.943 

2.2xl0 6 

CFD [19] 

0.537 

0.054 

0.943 

2.2xl0 6 

CFD [20] 

N/A 

0.056 

0.934 

1.8xl0 6 
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(a) Experiment (Re D =3.8M) [17] 


(b) CFD (Re D =2.2M) [19] 


Pt: 0.84 0.85 0.86 0.87 0.88 0.89 0.9 0.91 0.92 0.93 0.94 0.95 0.96 0.97 0.98 0.99 



(c) Present CFD-1 (Re D =3.8M) (d) Present CFD-2 (Re D =3.8M) 

Fig.A5 Comparison of total pressure contours at AIP (M^O.85) 
(See Fig. A3 for definitions of CFD-1 and CFD-2) 
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